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ABSTRACT 

The  acoustic  noise  spectrum  associated  with  traveling 
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Observed  noise  spectra  showed  that  low  frequency  noise 
can  be  modeled  by  incompressible  theory.  Hjqh  frequency 
noise,  apparently  resultinq  from  a  shock  wave*  can  be 
modeled  by  compressible  theory.  The  spectrum  was  seen  to 
shift  toward  lower  frequencies  than  predicted  possibly  due 
to  asymmetric  bubble  collapse. 

The  spectral  energy  per  bubble  was  experimentally 
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The  collapse  peak  pressure  amplitude  distribution*  the 
maximum  bubble  radius  distribution  and  the  nuclei 
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Chapter  l 


INTRODUCTION 

l.I  Historical  Background 

Cavitation  is  the  rupture  of  a  liquid  or  a  liquid-solid 
interface  caused  by  reduction  of  local  static  pressure  below 
the  vapor  pressure  of  the  fluid  medium*  This  rupture 
manifests  itself  as  a  macroscopic  bubble*  Microscopic  voids 
present  in  liquids  act  as  seeds  or  nuclei  for  cavitation 
bubbles  ( 1 } • 

Cavitation  is  a  prominent  source  of  noise  in  underwater 
acoustics.  It  is  basically  a  monopole  source  due  to  volume 
changes  but  some  dipole  and  quadrupole  components  are 
present  from  asymmetries  in  the  bubble  shape* 

When  bubbles  created  by  hydrodynamic  cavitation  are 
moved  by  the  flow  out  of  the  reqion  where  the  static 
pressure  is  less  than  vapbr  pressure*  the  vapor  condenses 
quickly  as  the  pressure  rises  resultinq  in  the  collapse  of 
the  bubble*  The  violence  of  this  collapse  produces  a  hiqh 
amplitude  pressure  wave  which  propagates  into  the  medium  as 
sound*  The  collective  intensity  of  thousands  of  bubble 
collapses  per  second  for  well  developed  traveling  bubble 
cavitation  is  such  that  a  major  noise  source  is  produced* 

The  study  of  the  mechanism  of  noise  production  by  a 


and 


collapsing  bubble  has  resulted  in  several  theories 
models  based  upon  numerous  assumptions  since  Besant  -(2) 
first  discussed  the  calculation  of  the  pressure  field  caused 
by  the  collapse  of  a  spherical  cavity  in  an  incompressible 
fluid  in  1859. 

The  first  significant  breakthrough  came  in  1917  when 
Lord  Rayleigh  (3)  calculated  the  pressure  field  around  a 
vacuous  cavity  and  predicted  the  collapse  time.  But.  because 
the  cavity  was  assumed  to  be  a  void*  an  infinite  wall 
collapse  velocity  was  predicted  as  the  cavity  radius 
approached  zero. 

The  next  logical  step  was  to  introduce  the  effect  of 
vapor  inside  the  cavity  recognizing  it  as  a  source  of 
limitation  to  the  wall  acceleration.  as  the  wall  velocity 
does  approach  the  speed  of  sound  however*  the  effect  of 
liquid  compressibility  has  to  be  accounted  for.  In  1952* 
Gilmore  (4)  used  the  K i r kwood-Bethe  approximation  to  model 
the  final  collapse  period  where  liquid  compressibility  is  a 
factor.  Mellen  (5)  substituted  the  use  of  Gilmore’s  theory 
for  final  collapse  and  calculated  the  resultant  shock  wave 
pressure  as  a  time  function.  The  development  of  a  shock  wave 
theory  was  furthered  by  Willard  (7).  He  also  observed  that  a 
produced  shock  wave  would  influence  the  growth  of  nearby 
nuc lei. 

The  effect  of  surface  tension  on  the  collapse  was  shown 
by  Bahl  and  Ray  (6)  to  be  significant  only  when  the  maximum 
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bubble  radius  was  of  very  small  magnitude.  They  found  that 


I  this  effect  was  not  significant  in  water* 

I 

1  As  theories  developed  for  predicting  bubble  collapse 

!  times  and  velocities*  the  next  major  area  of  investigation 

was  the  sound  power  spectrum  produced*  Mellen  (8) 
experimentally  derived  a  spectrum  and  in  1956  Fitzpatrick 

I 

and  Strasberg  (9)  predicted  a  peak  in  the  noise  spectrum  at 

I 

|  the  frequency  determined  by  the  reciprocal  of  the  Rayleiqh 

collapse  time.  Also*  they  used  the  Rayleiqh  analysis  to 
derive  the  far  field  acoustic  pressure  as  a  function  of  the 
second  time  derivative  of  the  bubble  volume.  Blake*  et.  al • 
(10)  used  these  single  bubble  analyses  to  scale 
experimentally  derived  spectra. 

j  More  recent  investigations  have  focused  on  the  role  of 

j  the  nuclei  in  cavitation  inception  and  scaling.  Van  der 

Walle  (ll)»  Apfel  (12)*  and  Lauterborn  (13)  have  all  tried 
to  explain  cavitation  inception  seal inq  in  terms  of  nuclei 
composition  and  dynamics.  Keller  (14)  and  later  Gates*  et. 
al.  (15)  used  1 iqht  scattering  and  holographic  methods  to 
study  the  relationship  of  nuclei  distribution  to  cavitation 
inception.  Gldenziel  (16)  did  a  study  of  what  he  called 
cavitation  susceptibility  and  of  cavitation  bubble  size. 

Models  qrew  more  complex  as  i deal i zat i ons  were  replaced 
by  natural  parameters.  Numerical  solutions  became  possible 
as  computers  were  made  more  powerful  and  affordable. 

|  Hickling  and  Plesset  (17)  did  a  numerical  analysis  of  the 


collapse  of  a  gas  filled  bubble  in  a  compressible  fluid  in 
1964*  Morozov  (18)  used  a  statistical  analysis  to  determine 
the  conditions  under  which  sinqle  bubble  theory  could  be 
used* 

The  most  recent  adaptation  of  sinqle  bubble  analysis  is 
the  investigation  of  the  effect  of  the  collapse  of 
asymmetric  bubbles*  Since  cavitation  often  occurs  near  a 
solid  boundary*  this  is  of  s i gn  i  f  i  cance*  Rayleiqh  assumed  a 
spherical  cavity  and  axisymmetric  collapses*  But  bubbles 
collapse  in  a  toroidal  shape  forming  water  jets  when  in  the 
proximity  of  a  solid  surface*  This  phenomenon  has  been 
investigated  by  Plesset  and  Chapman  (19),  Hsieh  (20), 
Kozyrev  (21),  and  Mitchell  and  Hammitt  (22).  The  use  of  hiqh 
speed  computers  has  resulted  in  models  which  aqree  well  with 
bubble  dynamics  recorded  by  high  speed  film.  Chahine*  et. 
al*  (23)  incorporated  asymmetric  bubble  collapse  into  a 
correlation  of  noise  and  cavitation  bubble  dynamics  in  1979. 

An  investigation  of  bubble  dynamics  and  the  use  of 
single  bubole  theory  to  experimentally  determine  the  noise 
per  bubble  collapse  was  conducted  by  Hamilton  (24).  He  found 
experimental  data  which  indicated  that  a  combination  of 
theories  was  necessary  to  describe  the  bubble  dynamics  and 
sound  power  spectrum* 

1*2  Objectives  of  Study 

Most  previous  cavitation  studies  have  focused  on  either 
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single  bubbles  generated  by  sparks  or  developed  body 
cavitation  types  such  as  vortex  and  sheet*  This  study  i-s  a 
blend  of  these  approaches.  Traveling  bubble  cavitation  lends 
itself  to  modeling  by  single  bubble  theory  yet  also  is 
subjected  to  near-solid  collapse  effects  which  are 
prominent  in  many  natural  applications* 

The  noise  for  traveling  bubble  cavitation  has  been 
found  to  be  more  intensive  than  that  from  frothy  cavitation. 
The  intensity  of  the  noise  from  traveling  bubble  cavitation 
is  determined  by  the  collective  intensity  of  many  unique 
bubble  collapses.  8ecause  of  this,  sinqle  bubble  analysis  is 
justified.  However,  several  theoretical  models  exist  for 
predicting  the  noise  produced  by  the  collapse  of  a  single 
bubble.  This  study  will  investigate  the  appl i cabi 1 i ty  of 
each  model  over  the  frequency  spectrum  of  the  noise.  The 
role  of  asymmetric  collapse  of  the  bubble  due  to  proximity 
to  a  solid  boundary  will  also  be  discussed. 

This  study  comprises  five  areas  of  investigation. 

The  dependence  of  the  noise  upon  flow  variables  such 
as  free-stream  pressure  and  velocity  is  examined.  The  shape 
of  the  sound  power  spectra  is  examined  to  determine  the 
frequency  ranges  of  noise  production.  A  third  area  of 
investigation  is  a  determination  of  the  noise  contribution 
of  each  bubble  by  measuring  the  total  noise  level  and  the 
number  of  bubbles  involved.  Previous  investigations 


(3*32*35)  have  concluded  that  the  maximum  radius  of  the 
vaporous  cavitation  bubble  is  a  key  parameter  in  -the 
determination  of  the  spectral  density  of  the  emitted  noise* 
Therefore*  the  relationship  of  the  peak  pressure  level 
produced  by  a  collapsing  bubble  and  its  maximum  radius  is 
i nves t i gated*  Also*  several  recent  studies  have  shown  an 
important  relationship  between  the  nuclei  size  and 
cavitation  (11*12*14*15)*  Hence*  it  is  important  that  the 
fifth  area  of  study  be  the  examination  of  the  relationship 
between  the  experimentally  determined  maximum  bubble  radius 
distribution  for  a  given  flow  condition  and  the  nuclei 
distr ibut ion. 

The  overall  objective  of  this  study  is  to  apply  the 
idealized  theories  developed  from  studies  of  the  collapse 
and  emission  of  sound  by  sperical  bubbles  created 
artificially  by  spark  gaps  and  lasers  to  the  asymmetric 
bubbles  of  body  cavitation  encountered  in  practice* 
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Chapter  2 


THEORY 

2*1  Single  Bubble  Analysis 

The  analytical  techniques  used  to  date  to  study 
traveling  bubble  cavitation  are  predicated  upon  the  veracity 
of  scaling  single  bubble  dynamics  for  high  collapse  rates. 
In  order  to  justify  the  use  of  single  bubble  analysis*  the 
condition  of  independence  of  collapse  rates  must  be 
satisfied.  Therefore*  bubbles  should  not  exist  in  such  a 
proximity  that  they  have  an  effect  upon  each  other. 

The  statistical  treatment  of  cavitation  noise  has  been 
particularly  popular  in  the  Soviet  Union,  levkovskii  (25) 
found  that  the  distribution  of  time  intervals  between 
successive  collapses  was  exponential.  Given  an  exponential 
distribution  for  time  intervals  with  a  parameter  of  \  for 
example*  the  distribution  of  event  occurrences  will  be  a 
Poisson  distribution  with  mean  rate  X  (26).  Akulichev  and 
01*shevskii  (27)  studied  the  statistical  characteristics  of 
cavitation  noise  using  the  theory  of  Poisson  processes  to 
theoretically  develop  a  continuous  energy  spectra.  Also* 
based  upon  a  Poisson  distribution  of  random  collapses* 
Morozov  (18)  modeled  cavitation  noise  as  an  aperiodic  train 
of  pulses  of  varying  amplitude  and  occurrence  rate. 


19 


Therefore*  the  spectrum  of  a  sinqle  bubble  collapse  would  be 
similar  to  the  spectrum  resultinq  from  a  collection-  of 
random  collapses* 

Supported  by  experimental  data*  these  theories  are 
employed  to  justify  the  investigation  of  sinqle  bubble 
dynamics  and  noise  production  which  is  the  basis  of  this 
study* 

The  energy  density  spectrum  of  a  single  bubble  collapse 
can  be  derived  from  the  power  spectrum  of  the  cavitation 


noise*  Rice  (28)  showed  that  the  enerqy  density  spectrum  is 
proportional  to  the  power  spectrum  for  shot  effect  noise 
such  as  travel inq  bubble  cavitation.  He  demonstrated  that 
the  total  noise  power  is  equal  to  the  number  of  bubble 
collapses  per  second  times  the  energy  per  bubble* 

2*1*1  Incompressible  Theory 

The  first  solution  to  the  problem  of  determining  the 
fluid  response  to  the  annihilation  of  a  spherical  void  in  a 
compressible  fluid  was  done  by  Lord  Rayleiqh  (3)*  His 
solution  was  based  upon  several  assumptions  includinq: 

-  an  incompressible  fluid* 

-  uniform  pressure  and  temperature  in  the 
cavity* 

-  the  cavity  and  surrounding  flow  field  are 
spherically  symmetrical* 

-  heat  conduction*  surface  tension*  and 
viscosity  factors  are  negligible* 


-  body  forces  are  zero* 

Let  us  first  trace  the  growth  of  such  a  cavity*  The 
velocity  potential  for  a  spherically  expanding  surface  is 
given  by 


<t>  = 


R2R 


(1) 


where  R  is  the  bubble  radius  at  a  qiven  time  and  R  =  dR/dt 
is  the  wall  velocity.  For  an  irrotational  motion  of  an 
incompressible  fluid*  the  equation  of  motion  is  generally 
given  by 

'  .2 


grad 


[i+w-if-o  <2> 

where  P(r,t)  is  the  instantaneous  pressure  and  p  is  the  mass 
density  of  the  fluid*  In  the  fluid*  the  following  applies 
for  an  expanding  sphere: 
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The  effects  of  surface  tension  and  viscosity  are 
negligible  compared  to  the  terms  in  Equation  (3)*  The 
pressure  drop  due  to  surface  tension  is  qiven  by 


2o 


AP  = 


st 


(4) 


where  is  °st  the  surface  tension*  The  viscous  term  is 

4v  • 

T  R  (5) 

where  v  is  the  kinematic  viscosity* 

At  the  interface  of  the  cavity*  r=R*  Equation  (3) 
becomes 
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(6) 
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acoustic  wavelength.  The  inner  and  outer  flows  can  be 
described  by  matched  asymptotic  expansions  for  the  boundary- 
value  problem.  As  described  by  Pierce  (29).  both  expansions 
represent  solutions  of  the  Helmholtz  equation.  The  outer 
expansion  terms  must  also  satisfy  the  Sommer f i el d  radiation 
condition.  The  inner  expansion  for  a  compact  source  is  given 
by 

3  2 

s  _  _  i  ?r  i 

(13) 


—  -  D  •  V  —  +  £  Q.. 

r  r  ,  , 


j,k  J  3Xj3xk 

where  the  three  terms  are  the  monopole.  dipole.  and 
quadrupole  cent r i ou t i ons  respectively.  The  outer  expansion 
is  expressed  as 


S(t-r/c)  _  v  .  D(t-r/c)  +  *  32  Qjk(t  r/c) 

r  r  r 


(19) 


j  ,k  3x^3xk 

It  is  the  acoustic  field  which  is  described  by  the 
compressible  theory  expression  of  Equation  (12). 

It  has  been  found  useful  to  express  the  pressure  in 
terms  of  two  non-dimensional  numbers.  The  coefficient  of 
pressure  is  used  in  hydrodynamics  to  normalize  the  pressure 
near  a  body.  It  is  given  by 


P-Pm 

Cp  =  r~r 

2  pUco 


(15) 


where  P  is  the  local  pressure  at  the  point  of  interest  and 
P  and  U  are  the  free-stream  pressure  and  velocity 
respectively.  The  cavitation  number  is  defined  as 

P  -P 
00  v 

a  =  r~T  <l6> 

7T  PU 
2 

where  P  is  the  vapor  pressure  of  the  fluid. 
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After  a  period  of  rapid  initial  growth*  the  bubble 
grows  at  a  fairly  constant  rate  as  indicated  by  Equation 


(8).  Thus*  Equation  (8)  reduces  to 

•  T2  nl/2 

*  ■  k  (pv 


p) 


(IT) 


when  P^  is  replaced  by  the  local  pressure  P  that  would  exist 
if  the  bubble  was  not  present  and  P(R)  is  replaced  by  its 


dominant  component*  P 


Substitution  of  non-dimensional 


terms  in  Equation  (17)  gives 


R  =  U 


-  3  (Cp  +  o) 


>r 


(18) 


Assuminq  relatively  constant  qrowth*  a  maximum  radius* 

R  »can  be  determined  from  the  product  of  growth  rate  and  the 
M 

time  of  growth  defined  as  the  period  in  which  the  bubble  is 

in  the  region  P  <  P  •  This  time  is  defined  as 

v 

d 


8  il  (l-  Cp)1/2 


(19) 


where  d  is  the  linear  path  of  the  bubble  through  the  P  <  P 

v 

reqion  and  Cp  is  the  averaqe  coefficient  of  pressure  in  that 
region*  The  product  of  Equations  (18)  and  (19)  is 


=  d 
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Baiter  (30)  modified  Equation  (20)  to  account  for  inertial 
effects  after  the  bubble  leaves  the  low  pressure  region. 
This  gives 
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where  Col  is  the  averaqe  pressure  encountered  by  the  bubble 
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as  it  enters  the  high  pressure  reqion  that  will  initiate 
col  1  apse* 

Rayleigh  derived  the  wall  velocity  during  collapse  from 
equating  the  sum  of  the  potential  and  kinetic  energies 


during  collapse  to  the  potential  enerqy  of  the  bubble  at 

maximum  volume*  This  yielded 

j  ttR3P  +  2npR2R3  =  j  ttrJp  .  (22) 

This  gives  a  wall  velocity  of 

r  fp3  v- r  1/2 

R  =  4  -  ^  -  1 

3  p  R3  (23) 

_  y_L 

For  other  than  the  initial  and  final  stages  of  collapse* 


Equation  (23)  simplifies  to 

V  rRiV/2  r 

•  2  P  M 

R  =  T  =  U  2 

3  p  R  00  . 
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Rayleigh  found  the  time  of  collapse  to  be 

Tc  -  °-915  V  Jy  *  °-915  £  ■fi’-  ^  •  «25' 

The  limit  of  Equation  (24)  as  the  bubble  collapses 
towards  a  zero  radius  is  seen  to  be  infinite  velocity* 
Physical  reality  dictates  otherwise*  The  presence  of  a 
permanent  gas  in  the  bubble  results  in  use  of  the  wall 
kinetic  enerqy  to  compress  the  gas  resulting  in  a 


termination  of  collapse  before  R  =  0* 


Assuminq  adiabatic  compression  of  the  gas*  Noltinqk  and 
Neppiras  (31)  derived  an  equation  for  the  wall  velocity  of  a 
partially  qas  filled  bubble  in  an  ideal*  incompressible* 
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the  minimum  radius  were  used  to  derive  the  volume 
acceleration  during  the  bubble  collapse  phase*  That  yielded 
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Substitution  of  Equation  (33)  into  Equation  (12) 
i nstantaneous  acoustic  pressure  as 


P(t)  = 


3rR2 


1-4 


R 


3^ 


The  peak  positive  acoustic  pressure  is 

P 

max  3  rR  2 


rR 
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(33) 


gives  the 


(34) 


(35) 


Less  than  1%  of  the  radiated  energy  is  due  to  the 
bubble  growth  phase  (1);  hence*  Equation  (34)  is  considered 
to  be  a  reasonable  approx i mat i on  of  acoustic  pressure 
produced  by  spherical  cavitation  bubble  collapse  in  an 
incompressible  fluid* 

The  acoustic  spectral  energy  density  is  found  by  first 
substituting  Equation  (24)  into  Equation  (11)  to  give 

P(r,t)  = 


O-V'8(cp  +  o>0'V0-8 


P  ’  (36) 
followed  by  taking  the  Fourier  transform  and  squaring  its 
absolute  value  to  finally  yield 


S(f) 
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2. 1.2  Compressible  Theory 


Incompressible  theory  has  been  demonstrated  to-  be 
adequate  for  most  of  the  bubble  growth  and  collapse  time 
period.  But*  in  the  final  collapse  stage*  the  wall  velocity 
exceeds  the  speed  of  sound  so  the  liquid  can’t  be  considered 
incompressi Dl e  any  longer.  Also  the  presence  of  qas  filled 
nuclei  increases  the  effective  fluid  compressibility.  The 
effect  of  compressibility  is  to  limit  the  wall  motion  thus 
reducing  the  wall  velocity.  Compressibility  limits  the  fluid 
inertia  associated  with  the  collapse. 

Because  an  exact  analytic  solution  cannot  be  found* 
numerical  solutions  have  been  done.  Flynn  (32)  outlined 
several  approximations  which  modify  the  Rayl eiqh-Pl esset 
equation  for  the  effects  of  compressibility.  The  one  which 
applies  to  the  final  stage  of  collapse  is  the  Ki rkwood-Bethe 
approximation.  Using  this  approximation*  Gilmore  (4)  derived 
a  modified  Ray  1 e i gh-Pl esse t  equation  as  follows 


where  H  is  the  specific  enthalpy  and  C  is  the  speed  of  sound 
in  the  liquid  at  the  bubble  interface.  Flynn  (33)  reviewed 
the  numerical  solutions  done  by  Hicklinq  and  Plesset  (17) 
and  Ivany  and  Hammitt  and  found  that  the  Gilmore  equation 
was  a  good  approx i mat i on  to  experimental  observations.  The 
Hicklinq  and  Plesset  analysis  took  the  solution  beyond  the 
final  collapse  into  the  region  where  the  bubble  rebounds  and 


generates  a  shock  wave* 

The  wall  velocity  is  given  by 
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where  HM  is  the  specific  enthalpy  at  the  maximum  bubble 
radius  and  R  is  the  bubble  radius  when  the  collapse 
initiates*  Equations  for  the  specific  enthalpy  for  water  and 
the  speed  of  sound  at  the  bubble  wall  are  given  by  Cramer 


(34)  as 

P(R) 


(  dP  =  7_  (3001  atm) 


1/7 


{  (P (R)+3000  atm)6/7-(Poo+3000  atm)6/7}  (40) 


and 


C  =  (Cq  +  6H)1/2 


(41) 


where  CQ  is  given  for  water  as  1482  m/s. 


Esipov  and  Naugol*nykh  (35)  estimated  the  maximum  wall 
velocity  as 


R  =  -  0.6  C 
max 


2 

pc  J 
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(42) 


where  Pm  is  the  pressure  at  the  interface  at  the  minimum 
bubble  radius. 


Gilmore  derived  the  pressure  field  as 


P(r,t)  -  f  (P2  +  +  f 


j  PR2 


(43) 
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2  dP2 
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identifying  the  first  two  terms  as  the  nonradiated  induced 
pressure  and  the  last  term  as  the  radiated  acoustic 
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pressure* 


When  the  bubble  collapses*  inertial  forces  in  the 
surroundinq  fluid  develop  rapidly  increasinq  kinetic  enerqy 
that  is  either  stored  in  the  compressible  fluid  or  converted 
into  internal  enerqy  of  the  bubble  qas*  The  inward  wall 
motion  is  arrested  by  the  bubble  internal  pressure  and  part 
of  the  stored  enerqy  is  radiated  as  a  shock  wave  from  the 
rebounding  bubble  (36)* 

Benjamin  (37)  predicted  a  maximum  shock  wave  pressure 


of 
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(44) 

where  Ps  is  in  kilobars  (10“’  kPa)* 

Esipov  and  Naugol*nykh  (34)  developed  an  equation 
describing  the  shock  wave  emitted  considering 
compressibility  and  the  presence  of  a  noncondensable  qas  in 
the  bubble*  The  shock  wave  is  assumed  to  be  of  the 

-t/e 

exponential  form  Ps(r*t)  =  Ps  e  •  Assuminq  y=  4/3  for  water 
this  approx imation  is  given  by 
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The  minimum  radius  is  related  to  the  maximum  radius 


R 

m 


by 


(48) 


and 


the 


acoustic 


Mach  number 


is  qiven  by 
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These  equations  are  valid  in  water  for  Ps  >>  p  c  /  7  and 

3/4 

r  >>  Rm  exp  (  M  '  )  • 

The  generation  of  the  shock  wave  is  modeled  as  a 
pressure  impulse*  The  pressure  qiven  by  Equation  (45)  when 
squared  and  the  Fourier  transform  taken  yields  an  enerqy 
density  spectrum  of 


Thi  s 
at  a 


spectrum  is 
rate  of  f~2  • 
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approximately  flat 


to  1/  2 it  9  then 


(50) 

rolls  off 


2*1*3  Asymmetric  bubble  collapse 


So  far*  all  of  the  described  theories  have  assumed  a 
spherical  shape  for  the  collapsinq  bubble*  In  practice* 
interest  is  often  focused  on  cavitation  near  a  solid 
surface*  In  this  case*  the  effect  of  proximity  to  a  solid  is 
enouqh  to  cause  asymmetric  collapse  of  the  cavity*  Many 
investigations  have  been  concerned  with  the  water  jet  formed 
in  such  an  asymmetric  collapse  (21*38*39)*  It  is  believed 
that  this  is  a  primary  mechanism  of  cavitation  damage*  Very 


31 


little  has  been  done  to  study  the  effect  of  the  asymmetric 
bubble  dynamics  upon  the  production  of  cavitation  noise. -The 
neglect  of  the  study  of  nonspherical  bubbles  is  not  due  to 
lack  of  recognition  of  the  problem;  rathert  the  difficulty 
in  medeling  a  nonspherical  motion  has  prevented  the 
derivation  of  an  exact  solution.  Therefore  most  of  the  work 
done  in  the  area  of  nonspherical  bubble  dynamics  has  been 
numerical  analysis. 

Plesset  (40)  has  found  that  theory  exaggerates  the 

♦  2 

importance  of  the  liguid  inertia  term  (3/2  R  term  in 
Eguation  6)  for  the  final  stage  of  the  collapse  of  a  bubble 
near  a  solid  surface. 

The  modeling  of  the  nonspherical  collapse  dynamics 
beqan  with  assuming  sliqht  deviation  from  sphericity.  The 
results  could  be  linearized  by  assuminq  the  amplitudes  of 
the  spherical  harmonics  were  small.  The  next  step  was  to 
introduce  nonlinearities  from  large  deviations. 

Blake  and  Gibson  (41)  modeled  the  early  staqes  of 
collapse  by  adapting  a  distribution  of  sources  and 
derivatives  (e.q.  doublets)  along  the  axis  of  symmetry 
inside  the  bubble.  This  technique  failed  to  adequately 
describe  the  final  stages  of  collapse  because  of  a 
restriction  in  the  resolution  of  the  surface  shape  due  to  a 
decrease  in  the  center  line  length. 

Mitchell  and  Hammitt  (22)  and  Plesset  and  Chapman  (38) 
used  the  Marker  and  Cell  finite  difference  approach  to  model 
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this  final  collapse  stage*  Lauterborn  and  Bolle  (42)  found 
good  experimental  agreement  with  the  results  of  Plesset  -and 
Chapman  in  addition  to  the  formation  of  a  stronq  shock  wave 
on  the  collapse  of  a  bubble  in  proximity  to  a  solid 
boundary. 

Extremely  complex  theoretical  derivations  of 
approxi mati ons  for  asymmetric  bubble  collapse  dynamics  have 
been  done  by  Shima  (43)  and  Hsieh  (20).  Hsieh  obtained  an 
equation  governinq  the  growth  and  collapse  of  a  confocal 
ellipsoidal  bubble.  Shima  found  that  the  time  of  collapse 
increases  as  the  ratio  of  to  the  distance  of  the  bubble 
centroid  from  the  solid  boundary.  This  conclusion  was  also 
reached  by  Plesset  and  Chapman  (38)  as  they  found  the 
collapse  time  to  be  described  by 


Tc  =  *M 


AP 


1/2 


(51) 


Chahine  (23)  found  that  this  could  also  be  described  by 

.1/2 

0.915  Rcmax 

AP 


(52) 


where  Rcmax  is  half  of  the  lenqth  of  the  major  diameter  of 
an  ellipsoidal  bubble. 

Of  particular  interest  concerning  asymmetric  collapse 
is  the  distribution  of  energy  between  production  of  a  water 
jet  and  noise  production.  Mitchell  and  Hammitt  (22)  found 
that  half  of  the  energy  from  the  impulse  produced  by  the 
asymmetric  collapse  went  to  the  jet  formation  and  the  other 
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half  went  to  noise  effects  and  viscous  dissipation*  They 
also  observed  that  the  bubbles*  originally  ellipsoids  with 
the  major  axis  parallel  to  the  solid  boundary*  progressively 
chanqed  shape  through  a  spherical  shape  until  just  at 
collapse  initiation  when  they  were  ellipsoids  with  the  major 
axis  normal  to  the  solid  boundary. 

In  one  of  the  few  instances  of  correlation  between  the 
produced  noise  spectra  and  dynamics  of  an  asymmetric  bubble 
collapse*  Chahine  (23)  found  that  the  noise  spectrum  shifts 
toward  lower  frequencies  as  the  ratio  of  R  to  the  centroid 

M 

distance  from  the  solid  wall  increases* 

The  conclusion  to  be  drawn  from  this  discussion  is  that 
the  complexity  of  model inq  nonspherical  bubble  collapse  has 
inhibited  its  inclusion  in  noise  production  theory.  For  the 
purposes  of  this  study*  the  spherical  collapse  assumption 
will  be  accepted. 

2.2  Scaling  Factors 

One  purpose  of  model  testinq  is  to  find  functional 
relationships  for  scalinq  model  c harac ter i st i c s  to  full  size 
conditions.  These  relationships  involve  powers  of  fluid 
static  and  bubble  dynamic  parameters.  Experimental  data  is 
used  to  demonstrate  the  appl i cabi  1  i  ty  of  theoretically 
derived  scalinq  laws.  The  practicality  of  a  particular 
theoretical  model  can  then  be  evaluated. 

The  parameters  of  qreatest  importance  in  model  tests  in 
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a  water  tunnel  are  freestream  velocity  and  pressure*  The 
pressure  is  commonly  expressed  in  terms  of  cavitation  number 
and  coefficient  of  pressure.  Of  course*  the  dynamic  pressure 
denominator  of  both  of  these  nond i mens i onal  numbers  is  a 

function  of  U  • 

00 

The  soectral  energy  density  predicted  by  incompressible 
theory  was  given  by  Equation  (35). 


SCO  ~  H.2- V*6<Cr 


+  o)1'V0-4 


(53) 


By  setting  all  but  one  parameter  constant*  if  possible*  one 
can  find  the  exponential  dependence  of  the  spectral  energy 
upon  that  parameter.  In  this  case*  R^  is  also  a  function  of 
tr  which  makes  a  direct  relationship  between  S(f)  and 
difficult  to  find. 

The  maximum  in  the  spectrum  produced  by  a  bubble 
collapse  in  an  incompressible  fluid  was  predicted  to  occur 
at  f  =  1/  Tc.  Recalling  Equation  (23)»  this  peak  frequency 
is  found  to  be 


f  =  0.77  —  (C  +  o) 
p 


1/2 


(54) 


Compressible  theory  yielded  a  relation  for  spectral 
enerqy  density  of 

.2 


S(f)  = 


(P0) 


(55) 


1  +  (2tt6  f )" 


where  Ps  and  0  are  calculated  by  Equations  (44)  and  (45) 
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respectively.  For  values  of  P  <  1CT3  atmospheres*  the  m3c h 
number  term  becomes  i ns i qn i f i c ant  in  Equation  (45).  The  role 
of  Rm  was  minimized  due  to  its  logarithmic  contribution.  For 
hiqh  frequencies  (f  >  1/(2  ^3))»  Equation  (50)  becomes 

2  0  5-2 

S(f)  ~  RmPg  f  (56) 

o 

which  indicates  a  high  frequency  roll  off  as  1/f  .  For  this 

study*  this  roll  off  is  above  the  frequency  ranqe  of 
interest.  Therefore*  the  approx i mat i on  of  Equation  (50)* 
after  substitution  of  the  results  of  Equations  (46-4c»),  for 
lower  frequencies 

S(f)  -  R^Pg  (57) 

is  the  important  result.  Note  that  this  is  independent  of 
frequency  and  free-stream  velocity  if  <r  is  held  constant  in 
determining  RM  •  A  relationship  for  PG  still  needs  to  be 
found. 

Blake  (10)  experimentally  developed  a  functional 
relationship  for  spectral  enerqy  density  of 

S(f)  ~  tc  P(r ,t) 2  (58) 

_ 2 

where  P(r»t)  is  the  time  averaged  mean  square  of  the  sound 
level  from  the  collapse. 

The  predominant  parameter  in  all  of  these  functional 
relationships  is  Therefore  it  is  necessary  to  develop  a 

better  understanding  of  the  factors  which  determine  the 
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maximum  bubble  radius  such  as  nuclei  size  distribution  and 
asymmetric  collapse  of  the  bubble. 


Chapter  3 


PROCEDURE 

3.1  The  12  Inch  Diameter  Water  Tunnel 

This  study  of  the  bubble  dynamics  and  resulting  noise 
from  traveling  bubble  cavitation  was  done  at  the  Garfield 
Thomas  Water  Tunnel  of  the  Applied  Research  Laboratory  at 
The  Pennsylvania  State  University  (ARL/PSU).  The  primary 
test  facility  was  a  rec i rcul at i nq  water  tunnel  with  a  12 
inch  (0.305  m)  diameter  circular  test  section.  The  water 
velocity  is  continuously  variable  to  80  feet  per  second 
( fps) (24.4  m/s).  The  pressure  in  the  test  section  can  also 
be  continuously  controlled  from  3  pounds  per  square  inch 
absolute  (psia}(20.7  kPa)  to  over  60  psia  (413.7  kPa).  The 
temperature  of  the  water  can  be  varied  from  60*  to  100*  F. 
(15#-38°C»).  The  uniformity  of  the  flow  is  maintained  by 
honeycombs  and  filtration  through  25  micron  filters.  A 
schematic  of  this  facility  is  given  in  Figure  1. 

Clear  plexiqlass  windows  were  mounted  on  the  top  and 
two  sides  of  the  test  section.  The  outer  faces  were  flat 
while  the  inner  faces  were  curved  to  maintain  the  circular 
cross  section.  The  specific  acoustic  impedance  of  the 
plexiglass  closely  matches  that  of  water  which  results  in 
minimum  sound  reflection  and  maximum  transmission  at  a  water 
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Illustration  of  the  12  Inch  Diameter  Water  Tunnel. 


and  plexiglass  interface*  Total  reflection  is  assumed  at  a 
water  and  air  or  plexiqlass  and  air  interface*  In  order  .  to 
reduce  the  optical  distortion  when  a  laser  or  photographic 
equipment  was  used*  one  window  was  manufactured  with  a  flat 
inner  face.  This  window  was  desiqned  so  that  cavitation 
would  not  be  induced  due  to  a  disturbance  in  the  flow  along 
the  wall  of  the  circular  test  section* 

A  major  objective  of  this  study  was  to  investigate 
noise  character i sti c s  at  frequencies  below  10  kHz*  So  that 
standinq  waves  were  not  produced  at  desired  frequencies*  it 
was  necessary  to  have  the  receiving  hydrophone  qreater  than 
one  wavelength  from  the  source  of  noise  which  was  bubble 
collapse  on  the  model  body*  Thus*  to  enable  measurement  down 
to  2*5  kHz*  the  hydrophone  had  to  be  situated  approximately 
24  inches  (0.59  m)  away*  Since  the  test  section  has  only  a 
six  inch  (0.15  m)  radius*  a  bottomless  aluminum  tank  was 
attached  to  the  top  plexiglass  window  and  filled  with  water 
to  provide  a  hydrophone  well* 

3*2  Model 

The  model  used  in  this  study  is  a  Schiebe  body*  The 
contour  of  this  half  body  is  the  result  of  a  disk  source  in 
a  uniform  stream*  This  particular  body  was  chosen  because  it 
does  not  encounter  laminar  separation  and  it  is  most 
susceptible  to  traveling  bubble  cavitation  (44). 

The  body  has  a  maximum  diameter  of  2*0  inches  (51  mm) 


and  a  face  diameter  of  1.5  inches  (38  mm).  The  nose  is 
attached  to  a  tapered  afterbody  which  was  strut  mounted-  in 
the  center  of  the  circular  test  section.  This  configuration 
is  shown  in  Figure  2.  The  actual  model  used  was  the  one 
which  Hamilton  used  for  the  majority  of  his  tests.  Wax  had 
been  used  to  fill  mounting  holes  where  transducers 
originally  had  been  installed.  The  transducers  were  removed 
because  they  were  in  the  near  field  of  the  noise  source. 

The  pressure  distribution  alonq  the  surface  of  the 
model  was  measured  by  Hamilton  (24).  The  local  Pressure  is 
expressed  by  the  dimensionless  quantity  known  as  the 
pressure  coefficient.  This  is  defined  as: 


where  P  is  the  local  pressure  for  a  specific  point  on  the 
model*  Pa.  and  Uoo  are  the  Free  stream  static  pressure  and 
velocity  and  p  is  the  mass  density  of  the  fluid  medium. 

A  least  squares  polynomial  curve  fit  was  used  to 
qenerate  the  continuous  curve  presented  in  Fiqure  3.  The 
normalized  length*  s/a*  can  be  correctly  interpreted  as 
inches  of  arc  lenqth  alonq  the  model  surface  since  the  body 
maximum  radius  is  one  inch.  The  equations  defining  the 
pressure  coefficient  for  the  Schiebe  body  used  in  this  study 
are  given  in  Table  1. 
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Figure  3.  Interpolated  Values  of  Coefficient  of  Pressure 
versus  Surface  Position  on  Schiebe  Body. 


43 


TABLE  1 


EQUATIONS  DEFINING  C  ALONG 

P 

SURFACE  OF  SCIEBE  BODY 

i)  0  <  -  <  0.625 
—  a  — 

Cp  =  -69.621  (j)7  -  79.022  (|)6  +  177.71  (j)5 

-90.038  (-)4  +  13.037  (~)3  +  0.085253  (~)2 
a  a  a 

-  0.25309  (-)  +  1.0068 
a 

ii)  0.625  <  -  <  0.650 

—  a  — 

C  =  165.00  (-)3 
P  a' 

iii)  0.650  <-  <  0.745 

—a  — 

C  =  122.46  (-)2 
P  a 

iv)  0.745  <-  <  0.770 

—a  — 

C  =  717.29  (-)3 
P  a7 

v)  -  >  0.770 

a  — 

c  =  0.330 

p  - 

0.390  -  (-) 

a 


-  438.02  (~)2  +  342.68  (-)  -  83.389 

a  a 


-  176.80  (~)  +  62.786 

a 


-  1697.7  (~)2  +  1340.9  (-)  -  354.26 
a  a 
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3*3  Air  Content 


The  total  air  content  of  the  water  tunnel  was 
maintained  at  a  relatively  constant  value  of  10*0  parts  per 
million  (ppm)  on  a  molar  basis*  This  was  reqularly  measured 
usinq  a  Thomas  Van-Slyke  blood  gas  apparatus*  City  water* 
used  to  initially  fill  the  water  tunnel*has  an  air  content 
of  about  17  ppm*  By  lowerinq  the  pressure  in  the  tunnel 
below  atmospheric  pressure  and  usinq  as  low  a  flow  velocity 
as  possible  to  circulate  the  water*  free-s tream  air  bubbles 
and  some  other  bubbles  driven  out  of  solution  were  removed 
by  a  vacuum  pump  from  collecting  domes  on  the  tunnel.  This 
was  done  until  the  air  content  was  stabilized  in  the 
vicinity  of  10*0  ppm*  Since  many  engineering  applications 
involve  higher  air  contents*  this  investigation  was  done 
with  a  higher*  relatively  constant  content  of  10.5  ppm. 
Also*  the  pressure  ranqe  between  incipient  cavitation  and 
the  development  of  an  attached  cavity  is  directly 
proportional  to  the  air  content*  The  relatively  high  air 
content  maintained  in  this  investigation  allows  a  wider 
operatinq  range  for  travel inq  bubble  cavitation. 

3.4  Measurement  of  Pressure  and  Velocity 

The  pressure  and  velocity  in  the  test  section  were 
determined  from  the  static  and  total  pressures.  The  static 
pressure  was  obtained  using  a  Validyne  pressure  transducer 
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comparing  the  pressure  at  a  test 


section  wall  tap  to 


atmospheric  pressure. 

The  velocity  was  found  from  the  total  pressure  by  use 
of  Bernoulli's  equation  (Eq.  2  ) •  On  a  streamline 
P  *  1/2  p  1?  ■*  p  g  x2  =  constant 
Therefore  putting  a  stagnation  probe  in  the  same 
streamline  at  the  same  height  qives 

P  ♦  1/2  p  U2  s  P 

oo  oo  o 

or 

P  -  P  =  1/2  p  U2 

O  oo  oo 

A  CEC  pressure  transducer  was  used  to  compare  static 
pressure  to  stagnation  pressure. 

The  outputs  of  both  transducers  was  displayed  on  a 
Dymec  integrating  digital  voltmeter.  The  transducers  were 
calibrated  with  a  piston  device  in  a  one  square  inch 
cylinder  filled  with  water.  Weiqhts  were  centered  on  the 
circular  disk  attached  to  the  piston  and  spun  to  reduce  the 
effects  of  friction.  Various  weights  were  used  to  produce  a 
range  of  pressures. 


3.5  Measurement  of  CavitaMon  Inception  and  Desinence 


The  propensity  of  a  particular  set  of  flow  conditions 
to  produce  cavitation  is  measured  by  the  cavitation  number 
given  in  Equation  (16)  as: 


a 


(60) 
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where  P  and  U  are  the  free-stream  static  pressure  and 
velocity  respectively*  P^  is  the  vapor  pressure  and  p  is  the 
mass  density  of  the  liquids  For  the  purposes  of  this  study* 
the  temperature  of  the  water  was  assumed  to  be  constant 
thereby  fixing  single  values  for  p  and  .  This  leaves  only 
free-stream  static  pressure  and  velocity  as  variable 
parameters* 

The  test  matrix  consisted  of  running  a  ranqe  of 
pressures  for  each  of  seven  velocities.  The  velocity  started 
with  25*0  fps  (7.62  m/s)  and  increased  in  increments  of  2*5 
fps  (.76  m/s)  to  an  upper  limit  of  40.0  fps  (12.20  m/s). 

With  the  velocity  set*  the  pressure  was  slowly  lowered 

until  traveling  bubble  cavitation  was  first  seen  with  the 

aid  of  a  stroboscope  to  freeze  the  motion  of  the  bubbles. 

This  point  is  cavitation  inception.  The  static  pressure  is 

used  to  calculate  the  incipient  cavitation  number*  o^.  The 

pressure  was  progressively  reduced  until  an  attached  cavity 

formed.  The  cavitation  number  for  this  point  was  also 

calculated.  The  pressure  was  then  raised  until  the  traveling 

bubble  cavitation  disappeared.  The  pressure  at  this 

desinence  point  was  noted  for  later  computation  of  the 

desinent  cavitation  number*  O’  • 

d 


Because  of  the  presence 

of 

free- 

stream  bubbles 

at 

relatively  high  air  contents  and 

the 

erratic 

nature  of 

the 

Foxboro  pressure  control  system 

for 

the 

tunnel  * 

the 

determination  of  the  points  of 

i ncept i on 

and 

des i nence 

are 
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quite  subjective*  Hamilton  tried  other  optical  and 
acoustical  methods  of  measurinq  cavitation  but  concluded 
that  visual  measurinq*  despite  its  i naccurac i es *  was  still 
the  preferable  technique* 

3*6  Noise  Measurement 

Noise  measurements  were  made  with  a  Celesco  LC-10  lead 
zirconate  titanate  piezoelectric  hydrophone.  The  hydrophone 
was  suspended  in  a  water  filled  aluminum  tank  mounted  atop 
the  test  section*  \hereby  allowing  more  accurate 
measurements  at  lower  frequencies.  It  was  directed  toward 
the  center  of  the  Schiebe  body  apor ox i matel y  one  inch 
downstream  from  the  nose.  This  conf i gur at i on  is  depicted  in 
Fiqure  4. 

The  hydrophone  signal  was  amplified  and  high-pass 
filtered  at  2.5  kHz  through  two  cascaded  filters.  The  output 
was  then  monitored  by  one  or  more  devices  dependinq  on  the 
measurement  being  done. 

Spectral  analysis  of  cavitation  noise  was  done  by 
inputing  the  filtered  signal  to  a  Spectral  Dynamics  SD-360 
real  time  FFT  processor.  Fourier  transforms  were  done  to  the 
signal  in  ranges  of  2.5-10kHz*  2.5-20kHz*  and  2.5-100kHz. 
The  spectra  were  transferred  to  an  IBM  System  7  digital 
computer  and  subsequently  to  an  IBM  System  34  digital 
computer.  There*  a  hydrophone  sensitivity  spectrum  was 
applied  to  the  received  spectra  and  then  they  were  plotted 
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Figure  4.  Configuration  of  Hydrophone  in  Water  Tank 
Used  to  Measure  Cavitation  Noise. 
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on  a  Calcomp  model  748  flatbed  plotter.  Two  flow  velocities. 


30.0  and  37.5  fps»  were  investigated  over  a  ranqe  -  of 
cavitation  numbers.  The  schematic  diagram  of  this 
instrumentation  is  shown  in  Fiqure  5. 

The  sensitivity  of  the  receiving  hydrophone  was 
established  by  using  another  LC-10  as  a  sound  source  on  the 
surface  of  Schiebe  body.  The  calibrating  source  was 
positioned  approximately  where  cavitation  bubble  collapses 
would  normally  occur.  This  was  driven  with  3.0  Vrms  white 
noise.  Several  spectra  of  the  receiving  hydrophone’s  siqnal 
were  recorded  for  later  calibration  use. 

Reverberation  tests  were  also  conducted  with  the  LC-10 
source.  Tone  bursts  were  used  to  produce  a  response  of  the 
receiving  hydrophone  which  was  recorded  on  an  oscilloscope. 
Photographs  were  then  taken  of  these  responses. 

An  analysis  of  the  bubble  collapses  per  unit  time  and 
produced  noise  was  desired  to  derive  the  noise  per  bubble 
collapse  as  a  function  of  pressure.  The  filtered  hydrophone 
siqnal  was  put  in  the  SD-360.  displayed,  and  copied  for  time 
domain  analysis.  Thus,  collapse  pulses  could  be  recorded 
graphically  for  later  counting.  The  filtered  signal  was  also 
monitored  by  a  true  rms  dB  voltmeter  set  on  a  long 
integration  time  to  give  the  produced  noise  level. 
Measurements  were  made  over  a  range  of  cavitation  numbers 
from  traveling  bubble  inceotion  to  attached  cavity  inception 
for  both  30.0  and  37.5  fps.  The  static  pressure  was  read  at 
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Figure  5.  Schematic  Diagram  of  Instrumentation  for 
Cavitation  Noise  Spectra. 


the  same  time  that  the  time  samples  were  obtained*  Fiqure  6 
shows  the  data  acquisition  system  used* 

3*7  Recordinq  Bubble  Dynamics 

An  investigation  of  bubble  dynamics  of  traveling  bubble 
cavitation  and  the  role  of  ft  was  done  by  videotapinq  the 
cav i tat i on* 

A  Sony  AV-3400  camera  and  recording  system  were  used* 
The  Schiebe  body  was  backlighted  with  a  strobe  behind  a 
diffusing  screen*  The  best  imaqes  were  obtained  when  the 
bubbles  were  directly  below  the  model  against  a  light 
background.  Imaqes  were  also  taken  of  a  scale  against  the 
body  for  later  use  in  dimensioning.  The  equipment 
configuration  is  illustrated  in  Fiqure  7. 

One  hour  of  video  tape  was  taken  for  each  of  two  flow 
conditions;  a  flow  velocity  of  30.0  fps  and  cavitation 
numbers  of  <T  =  .80  and  cr  =  .68  set  the  conditions*  While  the 
static  pressure  remained  relatively  constant*  deviations 
were  called  out  for  recording  on  the  tape  soundtrack. 

3*8  Measurement  of  Gas  Nuclei  in  Water 

The  role  of  nuclei  density  and  size  in  cavitation 
inception  has  been  the  subject  of  many  studies  in  the  recent 
past  including  those  by  van  der  Walle  (II)*  Johnson  and 
Hsieh  (45)*  Gates  ( 4  €» )  *  Keller  (14)*  and  Gates* 
Bi 1 let*et*al • (  15) .  Before  the  development  of  holographic  and 


Figure  6.  Schematic  Diagram  of  Instrumentation  for 
Cavitation  Bubble  Count  and  Noise  Data. 
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Figure  7.  Schematic  Drawing  of  Video  Taping  Apparatus. 
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1 iqht  scatterinq  systems  for  nuclei  countinq*  inaccurate 
estimates  were  all  that  were  possible*  The  total  qas  content 
was  used  as  one  such  estimator*  but  one  could  not  know  the 
proportion  of  free  to  dissolved  air  let  alone  the  nuclei 
size  distribution*  Another  estimation  could  be  derived  from 
the  attenuation  of  an  acoustic  siqnal  across  the  test 
section.  A  high  nuclei  density  results  in  hiqh  attenuation 
due  to  the  impedance  mismatch  presented  by  the  bubble 
sur  face* 

One  method  used  to  determine  nuclei  distributions  at 
ARL/PSU  is  Fraunhofer  holoqraphy.  A  two  step  procedure  is 
necessary  to  develop  an  image.  First*  a  sample  volume  is 
illuminated  with  a  high  intensity*  collimated  beam  of 
coherent*  monochromatic  1 iqht  from  a  ruby  laser.  The 
interference  between  the  coherent  backqround  and  the 
par ti cl e-di f f r acted  radiation  exposes  photographic  film  in 
the  far-field  of  the  nuclei.  This  Fraunhofer  hologram  is 
chemically  developed  and*  in  the  second  step*  illuminated  by 
another  coherent  beam  such  as  that  used  for  exposure.  A 
three  dimensional  image  of  the  test  section  volume  is 
produced  which  can  be  used  to  count  and  dimension  nuclei.  A 
detailed  description  of  the  system  is  qiven  in  Appendix  A. 

This  technique  does  have  some  drawbacks  though.  Because 
of  the  sensitivity  of  the  laser  to  temperature*  the 
intensity  of  the  pulse  at  a  specific  settinq  would 
fluctuate.  Oue  to  the  requirement  of  total  darkness*  a 
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wooden  frame  had  to  be  constructed  around  the  water  tunnel 


and  covered  with  heavy  black  plastic  making  temoerature 
control  even  more  difficult.  This  enclosure  is  shown  in 
Fiqure  8.  Therefore*  several  holograms  needed  to  be  taken  to 
obtain  one  qood  enough  for  reconstruction.  Also*  the  time 
required  to  count  and  size  the  particles  is  considerable. 

But*  since  holography  is  currently  the  most  accurate 
method  available  for  obtaining  nuclei  distributions*  it  was 
used  to  obtain  holograms  for  the  two  flow  conditions  which 
were  video  taped.  The  holograms  were  taken  just  previous  to 
the  commencement  of  taping  for  each  condition. 


! 
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Chapter  4 


RESULTS 

4*1  Cavitation  Characteristics 

Classical  theory  predicts  the  cavitation  number  for 
incipient  and  desinent  vaporous  cavitation  to  be  the 
maqnitude  of  the  minimum  pressure  coefficient  of  the  model. 
From  the  definitions  of  pressure  coefficient  and  cavitation 
number,  one  finds  that  only  if  a  <  -Cpmin  will  there  exist  a 
region  on  the  model  where  the  local  pressure  is  less  than 
vapor  pressure  thus  allowing  the  growth  of  a  vaporous  bubble 
Therefore,  -Cpmin  is  used  as  a  first  approx i mat i on  of  the 
incipient  cavitation  number. 

But  cavitation  does  not  always  incept  at  this  simply 
predicted  value  because  of  secondary  scale  effects  such  as 
model  size,  material,  gas  content  of  the  water,  and  tunnel 
velocity.  For  slow  tunnel  speeds  and  hiqh  air  contents,  van 
der  Walle  (111  found  that  inception  of  qaseous  cavitation 
often  occurs  at  cavitation  numbers  higher  than  predicted. 
When  the  tunnel  velocity  is  increased  and/or  the  air  content 
is  lowered,  the  vaporous  cavitation  is  predominant  as  the 
static  pressure  is  lowered.  It  has  been  observed  to  incept 
at  cavitation  numbers  below  predicted. 

The  condition  for  disappearance  of  cavitation  as  the 
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static  pressure  (and  therefore  the  cavitation  number)  is 
increased  has  been  often  found  to  differ  from  that  -  of 
inception.  This  phenomenon  is  described  by  Hoi  1  and  Treaster 
(47)  as  cavitation  hysteresis. 

The  result  of  cavitation  inception  and  desinence 
measurement  is  shown  in  Figure  9.  Several  cavitation 
character istics  are  demonstrated  by  this  graph.  First*  with 
-Cpmin  for  the  studied  model  equal  to  1.03*  the  same  trends 
are  seen  which  van  der  Walle  described.  It  is  very  difficult 
to  distinguish  between  vaporous  and  gaseous  cavitation  when 
measuring  inception.  Secondly*  desinence  is  shown  to  occur 
at  higher  cavitation  numbers  than  inception.  The  degree  of 
hysteresis  or  difference  between  inception  and  desinence 
decreases  markedl y  with  increased  tunnel  velocity.  Also*  the 
condition  required  for  the  inception  of  an  attached  cavity 
is  shown  to  be  relatively  constant.  Finally*  char acter i st i c 
cavitation  numbers  are  proportional  to  the  gas  content. 

The  subjectivity  of  visual  cavitation  determi nat ion  is 
evidenced  by  the  maximum  standard  deviation.  As  mentioned 
before*  efforts  have  been  made  to  employ  cavitation  event 
counters.  However*  this  method  is  not  entirely  free  of  human 
judgment  as  the  detection  threshold  and  the  inception  event 
rate  must  be  chosen  by  someone.  Gates*  et.  al .  (15) 
encountered  different  types  of  cavitation  and  cavitation 
locations  and  subsequently  discontinued  use  of  a  Keller  (14) 
type  optical  event  counter  in  favor  of  visual  measurement. 
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Figure  9.  Cavitation  Numbers  for  Inception  and  Desinence 

of  Traveling  Bubble  Cavitation  and  Inception  of 

an  Attached  Cavity  as  a  Function  of  U  . 
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4 .2  Single  Bubble  Cavitation  Noise 

The  objective  of  the  noise  tests  was  to  obtain  the  far 
field  spectrum  due  to  a  cavitation  bubble  collapse  over  the 
2.5  -  100  kHz  frequency  range.  In  order  to  have  the 
hydrophone  at  least  one  wavelength  of  the  lowest  desired 
frequency  away  from  the  source*  the  separation  distance  of 
the  model  and  the  hydrophone  needed  to  be  24  inches  (.59  m). 
But*  to  obtain  reasonable  hydrophone  response  in  the  2.5  -  4 
kHz  ranqe*  the  hydrophone  had  to  be  mounted  12  inches  (.30 
m)  from  the  model. 

Reverberation  tests  showed  a  strong  signal  received  by 
direct  path  and  indirect  path  reflection  off  of  the  bottom 
of  the  test  section  and  several  weaker  reflections  off  the 
hydrophone  tank  walls.  The  time  delay  until  the  first 
response  pulse  was  received  was  equivalent  to  the  time 
necessary  for  a  sound  wave  to  travel  the  direct  path  (12”)» 
The  time  until  the  second  strong  pulse  was  equivalent  to  a 
wave  reflecting  off  the  bottom  of  the  test  section  and 
travel inq  to  the  hydrophone  (26M).  The  reverberant  field 
after  shuttinq  off  a  steady  source  showed  an  exponential 
decay  with  two  modulations  in  the  decay  envelope.  These 
became  important  in  later  interpretat ion  of  the  bubble 
collapse  pulses.  Figure  10  shows  these  responses. 

Bubble  collapse  noise  bursts  were  recorded  in  the  time 
domain  for  later  analysis.  A  typical  plot  of  a  string  of 
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Figure  10.  Hydrophone  Response  to  a.)  1  Cycle  of  a 
25  kHz  Sinusoid  and  b.)  16  Cycles  of  a 
25  kHz  Sinusoid  Source  Input  on  a  Model 
Surface  at  Approximate  Location  of  Bubble 
Collapse. 


collapse  pulses  is  given  in  Fjqure  11.  The  plots  show  the 
nulls  evidenced  by  the  reverberation  tests  as  wel>  as 
strings  of  closely  packed  collapses.  These  groups  of  noise 
bursts  make  recognition  of  individual  pulses  difficult.  This 
is  one  source  of  possible  error  in  obtaining  the  bubble 
collapse  rate  for  each  flow  condition.  Discrimination  of 
discrete  pulses  is  most  difficult  at  low  cavitation  numbers. 

The  just i f i ca t i on  for  modeling  traveling  bubble 
cavitation  as  a  series  of  discrete  bubble  collapses  is 
dependent  upon  the  independence  of  collapses.  As  discussed 
in  Section  2.1,  a  Poisson  distribution  of  collapse  events 
satisfies  this  independence  criterion.  The  existence  of  this 
Poisson  distribution  can  be  established  by  the  existence  of 
an  exponential  distribution  for  the  time  intervals  between 
collapses.  Figure  12  shows  the  distribution  of  collapse 
intervals  obtained  from  a  200  msec  plot.  It  is  indeed  an 
exponential  distribution. 

With  collapse  independence  established,  the  next  steps 
toward  deriving  a  noise  per  bubble  relation  to  cavitation 
number  were  the  recording  of  produced  noise  level  and  the 
countinq  of  the  cavitation  pulses  on  the  time  domain  records 
to  establish  collapse  occurrence  rates.  Figure  13  is  a  plot 
of  the  rate  of  collapses  versus  cavitation  number.  The 
number  of  collapses  is  normalized  by  the  body  radius  divided 

by  the  free-stream  velocity  (a/U  )  resulting  in  a  Strouhal 

00 

number.  Normally  this  number  is  associated  with  the 
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Figure  12.  Noise  Burst  Interval  Distribution 
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frequency  of  shed  vortices.  Here  it  is  indicative  of  the 
frequency  of  another  time  deDendent  flow  phenomenum 
cavitation  bubbles.  Notice  that  at  low  cavitation  numbers* 
the  number  of  bubble  collapses  per  second  appears  to  be 
leveling  off.  This  may  be  due  in  part  to  the  difficulty  in 
distinguishing  individual  collapse  pulses.  The  collapse 
events  are  merginq  as  the  transition  from  traveling  bubble 
cavitation  to  attached  rinq  cavitation  is  approached.  The 
overall  cavitation  noise  level  as  monitored  on  a  true  rms  dB 
voltmeter  is  plotted  against  cavitation  number  in  Figure  14. 

The  sound  pressure  was  normalised  by  the  dynamic  pressure 

2 

(1/2  p  )  and  plotted  as  dBrms  re  1  jjPa.  The  least  squares 
fitted  slopes  for  both  Figures  13  and  14  are  for  the  ranqe 
•  67  <  <r  <  .87*  which  covers  most  of  the  traveling  bubble 

cavitation  ranqe.  Note  also  that  both  Fiqures  13  and  14  show 
a  possible  Reynolds  number  dependence. 

Division  of  the  total  noise  level  by  the  number  of 
bubble  collapses  results  in  the  noise  per  bubble  relation 
depicted  in  Figure  15.  Because  the  data  in  Figures  13  and  14 
are  relatively  linear  on  the  log  scales*  an  exponential 
relation  of  the  form  <yk  was  assumed.  The  least  squares 
obtained  exponents  are  given  in  Table  2.  The  curve  in  Fiqure 
15  was  obtained  by  subtracting  Figure  13  from  Fiqure  14.  As 
can  be  seen*  these  exponential  values  are  seemingly 
independent  of  velocity.  Over  the  measured  ranqe  of  O'  ♦  the 
data  shows  a  o  relationship  for  average  noise  oer  bubble. 
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Figure  13.  Nondimens ional  Number  of  Bubble  Collapses  per 
Second  versus  Cavitation  Number  for  Velocities 
U  =  30.0  fps  and  U  =  37.5  fps. 
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Figure  14.  Nondimens ional  Sound  Pressure  Level  versus 
Cavitation  Number  for  Velocities  =  30.0 
and  U  =37.5  fps. 
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Figure  15.  Cavitation  Noise  Level  per  Bubble  versus 
Cavitation  Number  for  Velocities  of  U 
30.0  fps  and  U  =  37.5  fps. 


TABLE  2 


FUNCTIONAL  RELATIONSHIPS  FOR  THE  DEPENDENCE 
OF  N,  L,  L„  ON  e 

D 


Functional  form  cr 


k 


assumed . 


Uoo  (fr  '  k 


A.)  Number  of  Noise  Bursts  per 


Second;  N 

30.0 

-11.23 

37.5 

-11.04 

B. )  Noise  Level;  L 

30.0 

-13.88 

37.5 

-13.58 

C.)  Noise  Level  per  Bubble;  Ln 

D 

30.0 

-  2.65 

37.5 

-  2.54 
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The  spectrum  analysis  of  the  cavitation  noise  was  done 
with  a  Spectral  Dynamics  360  analyzer  and  two  IBM  computers. 
A  Fourier  transform  was  done  on  the  siqnal  received  from  the 


hydrophone  to  produce  the  raw  spectrum.  In  order  to  derive 
the  true  energy  density  spectrum*  the  frequency  response  of 
the  hydrophone  must  be  accounted  for  al onq  with  the 
elimination  of  the  reverberant  field  caused  by  the  proximity 
of  the  tunnel  walls.  The  procedure  used  to  develop  such  a 
calibrated  sensitivity  spectrum  is  outlined  in  Appendix  B. 
This  spectrum  was  subtracted  from  the  raw  spectrum  to  yield 
the  enerqy  density  spectrum.  Low  frequency  noise  was  removed 
from  the  spectra  by  high  pass  filtering  of  the  input  siqnal 
at  2.5  kHz.  The  filter  used  has  a  96  dB  per  octave 
attenuation  rate  below  the  high  pass  point*  while  above  this 
point  unity  gain  exists. 

The  resultant  spectrum  shows  many  resonances  which  may 
be  due  in  part  to  calibration  techniques  and  standing  waves 
in  the  water  tunnel.  The  cal ibration  set-up  placed  the 
source  hydrophone  on  the  top  of  the  model  at  the  approximate 
position  of  bubble  collapse.  However*  the  actual  collapses 
are  distributed  as  a  ring  source  around  the  model.  The 
sensitivity  spectrum  will  not  remove  characteristic 
resonances  from  such  a  distribution.  Also*  standing  waves 


are  present  in  the  water  tunnel  which  add  peaks  at  4860  Hz 
and  its  multiples.  This  same  phenomenon  has  been  encountered 
in  other  tests  at  this  facility.  The  peaks  are  removed  by  an 


averaginq  process*  This  technique  involves  drawing  an 
envelope  around  the  jaqged  spectrum  as  shown  in  ^iqure  16. 
The  upper  and  lower  values  of  the  envelope  are  converted  to 
pressure  values,  averaged,  and  changed  back  to  a  decibel 
value  as  the  averaged  maqnitude  for  that  frequency. 

Four  samples  were  taken  for  each  flow  condition  to 
determine  repeatability.  The  spectral  levels  were  all  within 
a  four  dB  band.  This  variation  is  due  primarily  to  slight 
differences  in  the  cavitation  numbers  because  of  the 
difficulty  in  exactly  repeating  a  specific  tunnel  pressure. 

Averaged  spectra  for  typical  cavitation  numbers  are 
shown  in  Figures  17  and  18  for  free-stream  velocities  of 
30.0  and  37.5  fps.  The  frequency  range  is  2.5  to  80  kHz 
since  the  spectrum  analyzer  has  an  aliasinq  filter  at  the 
80%  spectrum  point  to  prevent  interference  of  hiqh  frequency 
and  low  frequency  components.  The  low  pass  filter  reduces 
the  amplitudes  of  the  hiqh  frequencies  to  avoid  the  effects 
of  aliasing. 

For  higher  frequencies.  the  spectrum  is  relatively 
level  and  has  a  weak  dependence  upon  frequency.  The  spectrum 
seems  to  have  reached  a  plateau  at  the  higher  frequencies. 
The  higher  frequency  roll-off  point  is  predicted  by 
compressible  theory  but  it  is  beyond  the  measured  frequency 
range. 

The  absolute  sound  pressure  levels  indicated  in  Figures 
17  and  18  are  lower  than  expected  possibly  due  to 
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Averaged  Traveling  Bubble  Cavitation  80  kHz 
Noise  Spectra  for  U  =  30.0  fps. 


Figure  18.  Averaged  Traveling  Bubble  Cavitation  80  kHz 
Noise  Spectra  for  U  =37.5  fps. 


attenuation  by  bubbles  in  the  flow  near  the  top  of  the  test 
section.  This  condition  manifested  itself  mostly  at  low 
cavitation  numbers  and  especially  at  lower  free-stream 
vel oc i t i es. 


4.3  Low  Frequency  Characteristics 


Cavitation  noise  spectra  were  also  measured  over  2.5  - 
10  kHz  and  2.5  -  20  kHz  ranges  for  better  observation  of  low 
frequency  noise.  Figures  19  and  20  are  typical  20  kHz 
spectra  for  the  tested  velocities.  Once  again*  80%  of  the 
spectrum  is  plotted  because  of  the  aliasinq  filter.  These 
fiqures  indicate  a  minimum  difference  between  cavitation  and 
background  flow  noise  as  infered  by  the  <r  =  1.03  curves 
between  7  and  8  kHz.  3elow  this  frequency*  there  is  another 
hump  of  noise  enerqy. 

Greater  low  frequency  detail  is  given  in  Fiqures  21  and 
22.  Note  that  the  spectrum  below  2.5  kHz  has  been  filtered 
out  .  But*  it  appears  that  there  is  a  peak  in  low  frequency 
enerqy  at  a  frequency  between  5  and  6  kHz  which  corresponds 
to  the  reciprocal  of  the  collapse  time*  tc  • 
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eraged  Traveling  Bubble  Cavitation  16  kHz 


Figure  21.  Averaged  Traveling  Bubble  Cavitation  8  kHz 

Noise  Spectra  for  U  =30.0  fps. 
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Averaged  Traveling  Bubble  Cavitation 
Noise  Spectra  for  U  =37.5  fps. 


4.4  Analysis  of  Bubble  Dynamics 

Video  taping  was  used  to  record  a  number  of  cavitation 
bubbles  distributed  over  the  Schiebe  body.  Video  tapes  were 
recorded  for  two  conditions:  O'  =  .80  and  <r  =  .68  with 
free-stream  velocity  of  30.0  fps.  The  total  air  content  was 

8.5  ppm  on  a  molar  basis. 

The  bubbles  were  measured  along  the  backlighted  profile 
of  the  model.  Measurements  were  taken  of  the  depth  and  width 
of  the  vapor  bubbles  as  well  as  their  positions  on  the  body. 
Previous  films  of  traveling  bubble  cavitation  have  shown  the 
bubbles  to  be  elliotical  in  profile  yet  circular  when  viewed 
normal  to  the  model  surface.  Therefore  the  bubble  can  be 
considered  an  oblate  spheroid  for  the  purpose  of  calculating 
volume.  The  depth  is  the  minor  semiaxis  and  the  width  is  the 
major  axis.  The  bubbles  tended  to  be  flattened  out  at  hiqher 
flow  velocities. 

Using  the  calculated  volume*  the  radius  for  an 
equivalent  volume  sphere  was  found  because  the  cavitation 
collapse  is  being  modeled  as  a  spherical  collaose.  Recall 
from  Section  2.1.3  that  Mitchell  and  Hammitt  observed  that 
asymmetric  cavitation  bubbles  tended  to  alter  in  shape  from 
the  major  axis  parallel  to  the  model  surface  to  the  major 
axis  beinq  perpendicular  to  the  surface.  Analysis  of  the 
video  tapes  showed  a  similar  trend  of  the  flattened  bubbles 
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becominq  more  spherical  just  prior  to  collapse* 

Because  the  imaging  rate  of  a  video  tape  system  is 
relatively  slow  (  about  30  images  per  second)  only  one  image 
of  each  bubble  was  seen.  Therefore  a  larqe  sample  of  bubbles 
was  taken  py  measuring  size  and  location  of  each  as  seen  on 
the  video  monitor. 

The  ultimate  objective  was  to  find  a  distribution  of 
sizes  of  bubble  radii  at  the  maximum  radius  position  on  the 
model.  The  point  on  the  model  where  the  bubble  will  be  at 
its  maximum  volume  was  determined  by  developing  histograms 
of  bubble  sizes  at  several  points  alonq  the  model.  The  point 
associated  with  the  histogram  having  the  highest  mean  radius 
was  assumed  to  be  the  position  of  .  Further  analysis  of 
the  video  tapes  concentrated  on  the  .1  inch  region  around 
that  point. 

The  rate  of  occurrence  of  a  particular  bubble  size  was 
scaled  from  the  results  of  the  bubble  size  distribution  for 
the  narrow  surface  area  on  which  the  bubble  profile  could  be 
seen.  The  scaling  included  the  effect  of  the  difference  in 
static  pressure  around  the  body.  Lower  static  pressure  on 
the  top  of  the  model  allowed  bubbles  there  to  grow  larger 
than  they  would  have  on  the  bottom  of  the  model. 

The  resultant  maximum  bubble  radius  distributions  are 
qiven  in  Fiqures  2 3  and  24  •  The  averaqe  maximum  bubble 
radius  was  .388  mm  for  the  cr  =  .80  condition  with  standard 
deviation  .194  mm.  The  averaqe  maximum  bubble  radius  for  the 
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Figure  23.  Maximum  Cavitation  Bubble  Radius  Distribution 

for  U  =  30.0  fps  and  o  =  .80. 
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<T  =  .68  condition  was  .543  mm  with  standard  deviation 
mm.  The  sample  size  of  bubbles  at  the  point  of  maximum 
bubble  radius  is  smaller  for  the  <r  -  .80  condition  than  for 
the  lower  cavitation  number  condition.  This  was  due  to  the 
decreased  rate  of  cavitation  events. 

This  technique  for  studying  bubble  dynamics  is  limited 
in  that  it  is  a  statistical  approach.  it  cannot  be  used  to 
trace  the  growth  and  collapse  of  individual  bubbles  which 
would  be  necessary  to  investigate  asymmetric  collapse.  The 
poor  resolution  of  the  TV  screen  (♦  .09  mm)  also  limited  the 
accuracy  of  measurements. 
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4.5  Correlation  of  Bubble  Dynamics  to  Noise  Spectrum  Levels 

The  peak  pressure  amplitude  distribution  produced  by 
collapsing  bubbles  was  developed  from  the  time  domain 
records  produced  by  the  SD  360  analyzer.  The  dB  level  above 
background  was  measured  for  each  pulse  produced  by  a  bubble 
collapse.  A  histogram  was  used  to  collate  the  results. 
Finally.  a  curve  fit  was  used  to  provide  a  smooth 
distribution  curve.  Figures  25  and  26  show  the  distribution 
of  peak  pressure  amplitudes  produced  by  collapsing  traveling 
bubbles.  Note  that  the  distribution  is  loqnormal  as  was  the 
R  w  di str i bution. 

A  comparison  of  the  distribution  to  the  peak 

amplitude  distribution  yields  the  results  in  Figure  27. 

Equation  (35)  predicts  a  relationship  of  P  -v  R  for  bubbles 

M 

collapsing  far  from  the  surface.  The  relationship  is 

approximately  linear.  But*  it  cannot  be  known  how  many 

pressure  peaks  were  produced  by  bubble  rebounding  nor  if  a 

particular  magnitude  R  is  more  susceptible  to  rebound  than 

M 

another . 
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Peak  Pressure  Amplitude  'dBVrms  above  Background) 


Figure  26.  Cavitation  Bubble  Collapses  per  Second  versus 
Peak  Pressure  Amplitude  for  U  =  30.0  fps 
and  a  =  .68. 
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Figure  27.  Peak  Pressure  Amplitude  versus  Maximum 

Bubble  Radius  for  U  =  30.0  fps  and 
,  .  00  r 
CT  =  .68. 
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4.6  Correlation  of  Bubble  Dynamics  to  Nuclei  Distribution 


Nuclei  di s tr i but i ons  were  holographically  obtained  just 
prior  to  the  video  tapinq  of  each  of  the  two  tested 
cavitation  conditions.  The  distributions  are  for  a  volume 
just  upstream  of  the  Schiebe  body. 

The  results  of  the  analysis  of  the  holograms  are  shown 
in  Figures  28  and  29.  In  each  case*  the  distribution  appears 
to  be  lognormal.  Oldenziel  showed  the  form  of  this 
distribution  to  be 


nR(R)  = 


V 


-l/2s' 


/2tT  sR, 


exp 


0 


1  ,  2  R 

— 2  ln 
2s  R0 


(61) 


where  s  is  the  standard  deviation  and  n0  is  the  number  of 
bubbles  per  unit  volume. 

The  relationship  of  the  nuclei  distribution  to  the 
maximum  radii  distribution  is  not  altogether  clear.  It  may 
well  be  that  the  cavitation  vapor  bubbles  grow  from  the 
average  size  nuclei  while  the  small  nuclei  dissolve  from 
surface  tension  and  the  large  nuclei  become  stable 
gas-filled  bubbles.  One  observation  thouqh  is  that  a  greater 
percentage  of  nuclei  eventually  develop  into  cavitation 
bubbles  under  conditions  with  lower  cavitation  numbers. 
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TABLE  3 


BIN  SIZES  FOR  HOLOGRAPHIC  NUCLEI  DISTRIBUTION 


CHANNEL  BUBBLE  DIAMETER  (|»m) 


0 

2.0 

-  4.5 

1 

4.5 

-  7.5 

2 

7.5 

-  10.5 

3 

10.5 

-14.5 

4 

14.5 

-  18.7 

5 

18.7 

-24.0 

6 

24.0 

-  31.5 

7 

31.5 

-40.5 

8 

40.5 

-  52.0 

9 

52.0 

-  68.0 

10 

68.0 

-  86.5 

11 

86.5 

-  123 

12 

123 

-  195 

13 

195 

-  281 

14 

281 

-  467 

15 

467 

-  726 
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Chapter  5 
SUMMARY 


5.1  Discussion 

One  of  the  objectives  of  this  study  was  to  gauqe  the 
appl icabi 1 i ty  of  various  theories  for  predicting  noise 
produced  by  collapsing  bubbles  of  traveling  bubble 
cavitation*  This  has  taken  the  form  of  evaluation  of  the 
produced  spectrum  and  sound  pressure  level* 

As  seen  in  the  figures  of  Sections  4*2  and  4.3*  the 
cavitation  noise  spectrum  is  composed  of  two  portions  which 
are  modeled  by  different  theories.  The  low  frequency  portion 
is  described  by  incompressible  theory*  There  is  a  peak  in 
the  spectral  energy  density  at  an  approximate  frequency 
determined  by  the  reciprocal  of  the  collapse  time*  Tc.  The 
presence  of  a  solid  boundary  results  in  asymmetric  collapse 
of  the  bubbles  which  shifts  this  peak  to  1/(1* 5tc)*  The 
latter  frequency  was  observed  by  both  Blake  (10)  and  Ross 
(1)*  This  is  also  approximately  equivalent  to  Chahine's 
findinqs  given  by  Equation  (52)*  Blake  and  Gibson  (41)  also 
found  this  relationship*  Calculating  the  collapse  time  from 
Equation  (25)  and  using  1/(1. 5Tc )  to  find  the  peak*  the 
frequency  found  for  =  30.0  fps  and  <r  =  *68  was  4*9  kHz. 

Similarly  for  U  =  37.5  fps  and  <T  =  .68*  the  peak  frequency 
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derived  was  6.1  kHz.  When  these  values  are  compared  to  the 
peaks  in  Figures  19  and  20*  a  close  correspondence  can  he 
seen.  The  values  of  tc  derived  were  on  the  order  of  .1  msec 
for  the  Bayleiqh  equation  and  .14  msec  for  the  Gilmore 
equation  but  the  time  domain  records  show  a  Tc  of  closer  to 
•18  to  .20  msec  to  be  more  accurate.  It  is  reasonable  thouqh 
to  expect  1/Tc  to  overpredict  the  peak  frequency  because  its 
derivation  is  based  upon  the  assumption  of  a  void  cavity 
collapsing  symmetrically.  The  presence  of  a  permanent  gas  in 
the  bubble  will  retard  the  wall  motion  resulting  in  a  longer 
col  1  apse  time. 

This  low  frequency  spectral  shape  has  been  observed  by 
many  invest i qator s  including  Mellen  (8)»  Fitzpatrick  and 
Strasberq  (9)*  and  Barker  (48).  This  portion  of  the  spectrum 
extends  to  approx i mate  1 y  the  frequency  determined  by  the 
reciprocal  of  the  combined  qrowth  and  collapse  time  of  the 
bubble  reachinq  a  minimum. 

Above  that  point*  the  spectrum  rises  until  reachinq  a 
plateau  which  covers  most  of  the  10  -  100  kHz  range.  This 
can  be  seen  in  Figures  17  and  18.  This  level  spectrum  is 
indicative  of  a  shock  wave.  Compressible  theory  is  best 
suited  for  this  portion  of  the  spectrum.  Hamilton  (24)* 
Blake  *  and  Barker  all  found  a  similar  level  spectral  shape 
for  this  higher  frequency  range.  Lauterborn  and  Timm  (49) 
even  observed  two  shock  waves  produced  by  one  asymmetric 
collapse.  A  high  end  roll-off  point  is  predicted  by 
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compressible  theory*  but  the  measured  frequency  ranqes  do 
not  extend  that  high. 

Since  incompressible  theory  is  applicable  until  the 
final  staqes  of  collapse*  it  is  interesting  to  compare  the 
noise  level  per  bubble  results  of  this  study  to  the  levels 
predicted  by  incompressible  theory. 

Substitution  of  the  maximum  bubble  radius  relation  of 
Equation  (20)  into  the  incompressible  prediction  for 
spectral  energy  density  given  by  Equation  (53)  yields  the 
following  dependence  upon  a  : 

S  -  <<T  +  a)9/5(Cp'  +  a)6/5  (62) 

where  Cp*  is  the  average  coefficient  of  pressure  in  the  °<Pv 
reqion  and  Cp*  is  the  local  coefficient  of  pressure 
initiating  collapse.  As  an  estimate  of  Cp*  it  is  reasonable 
to  use  Cpmin  giving  Cp  =  -1.  Using  the  results  shown  in 
Figure  3  and  the  video  tape  results  for  point  of  collapse 
initiation*  a  value  of  Cp*  =  -.33  was  found.  When  these 
values  are  incorporated  into  Equation  (60)*  it  becomes 

S  ~  (1  -  a)9/5(a  -  .33)6/5  (63) 

This  is  plotted  in  Figure  30.  Also  plotted  are  the 
experimentally  derived  functional  relationships  found  by 
this  study  (  <T-^'^)  and  Hamilton  (  o1-®'^).  Note  that  each  is  a 
reasonable  approximation  to  the  curve  over  the  cavitation 
ranqes  tested.  This  suggests  that  incompressible  theory  can 


10  log  [(1  +  O  y,J  (  a  -  1/3) 
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Figure  30.  Comparison  of  Experimental  Results  of 

Hamilton  (24)  and  Present  Study  to  Incompressible 
Theory  for  Traveling  Bubble  Cavitation. 
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model  the  spectral  energy. 

The  absolute  sound  pressure  levels  are  lower  than 
expected  due  to  high  frequency  attenuation  by  bubbles  near 
the  top  of  the  test  section.  Another  contributing  factor 
suqgested  by  Chahine  (23)  is  the  directivity  of  noise 
produced  by  the  collapse  of  a  bubble  near  a  solid  boundary. 
He  found  that  noise  is  radiated  normal  to  the  surface  and 
attenuated  in  directions  parallel  to  the  surface.  This  is 
possibly  due  to  the  presence  of  many  small  bubbles  resulting 
from  the  bubble  collapse.  The  cylindrical  body  used  in  this 
study  would  tend  to  reflect  Chahine's  findings  because  of 
the  distribution  of  bubble  collapses  around  the  body.  Jones 
and  Edwards  (50)  found  a  peak  pressure  of  320  dB  re  1  pPa 
for  bubble  collapse.  This  is  about  10  dB  higher  than  the 
level  indicated  in  Figure  15. 

The  relationship  between  the  nuclei  distribution*  the 
maximum  bubble  radius  distribution*  and  the  peak  pressure 
amplitude  distribution  is  a  progressive  one.  All  of  these 
distributions  were  found  to  have  a  lognormal  shape.  The 
linear  relationship  between  the  peak  pressure  amplitude  and 
maximum  bubble  radius  contrasts  that  found  by  Chahine  of 

P  ~  .7  T2,5  (64) 

where  T  is  substituted  due  to  its  linear  dependence  on  R  . 

M 

It  should  be  realized  however  that  Chahine’s  results  were 
for  bubbles  collapsing  far  from  a  solid  boundary. 
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The  one  parameter  which  occurs  in  all  of  the  noise 
prediction  theories  that  is  the  most  difficult  to  determine 
is  the  maximum  bubble  radius*  Several  relations  have  been 
found  to  predict  RM*  Equation  (20)  predicts  for  =  30.0 
fps  and  or  =  *80  and  *68  to  be  *364  mm  and  .939  mm 
respectively.  31ake  et.  al.  found  the  relationship 


(-  o  -  C  .  ) 
2  pmin 


U  (1  -  C  ) 

oo '  p' 


for  travelinq  bubbles  on  a  solid  surface.  This  yields  .115mm 
and  .446  mm  for  the  respective  flow  conditions.  The  results 
of  video  tape  analysis  were  .388  mm  and  .543  mm  for  the  two 
cavitation  conditions  tested. 

The  lower  latter  two  sets  of  numbers  may  be  the  result 
of  the  influence  of  a  solid  boundary  on  bubble  growth. 
Therefore  it  is  necessary  to  investigate  the  dynamics  of 
qrowth  further  in  order  to  develop  a  more  accurate 
prediction  of  ^ . 


5.2  Conclusions 


This  study  has  examined  the  relation  of  bubble  dynamics 
to  the  noise  produced  by  traveling  bubble  cavitation.  The 
spectral  shape  and  sound  pressure  levels  have  been  used  to 
evaluate  the  ability  of  theories  derived  for  the  collapse  of 
single  bubbles  to  describe  well  developed  travelinq  bubble 


cav i tat i on. 


The  conclusions  to  be  drawn  from  the  study  can  be 


summarized  as  follows* 


*  For  frequencies  greater  than  approximate! y  10  kHz* 
the  spectral  enerqy  density  is  a  very  weak  function 
of  frequency  indicative  of  a  shock  wave* 
Compressible  theory  best  describes  this  portion  of 
the  spectrum* 

£  For  frequencies  below  10  kHz*  the  spectral  enerqy 
density  is  a  function  of  frequency  with  a  peak  at 
l/Tc •  This  is  indicative  of  the  noise  produced  by 
bubble  collapse  in  an  incompressible  liquid* 
Incompressible  theory  is  better  suited  for  this 
lower  frequency  range* 


*  The  magnitude  of  l/Tc  is  overpredicted  by  incom¬ 
pressible  theory  due  to  the  presence  of  a  permanent 
gas  in  the  cavitation  bubble. 

*  Asymmetric  bubble  collapse  may  be  responsible  for 
the  shift  of  the  spectrum  toward  lower  frequencies 
than  predicted  by  theories  derived  for  spherical 
bubble  collapse* 


*  The  exper imental 1 y  derived  dependence  of  noise 
level  per  bubble  collapse  of  l  ”  cr-2*6is  a  close 
approx i mati on  to  the  relation 


(C  + 


o)9/5(a 


v> 


6/5 


over  the  range  *67  <  a  <  *8A. 


$  The  amplitude  distribution  of  peak  pressure 

produced  by  a  collapsing  bubble*  the  maximum  bubble 
radius  distribution*  and  the  nuclei  distribution 
are  all  lognormal* 


This  study  has  touched  on  several  areas  which  require 


further  investigation*  The  effects  of  asymmetric  bubble 
collapse  upon  noise  production  is  important  to  most 
practical  applications  of  noise  scalina.  The  dependence  of 


the  maximum  bubble  radius  distribution  upon  the  free-stream 


velocity  and  cavitation  number  is  needed  since  RM  olays  such 
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an  important  role  in  theoretical  noise  predictions.  Also* 
the  hiqh  f requency  rol 1 -of f  is  another  area  to  be  studied. 
One  further  recommendation  for  followup  is  to  establish  the 
exponential  dependence  of  the  spectral  enerqy  density  upon 
the  cavitation  number  over  other  ranqes  of  <y.  with  the  peak 
pressure  amplitudes  from  the  time  domain  plots*  a  follow-on 
project  could  model  the  cavitation  noise  process  as  a 
Poisson  distribution  of  Dirac  delta  functions  of  variable 
amplitude  and  obtain  the  autocorr el ati on  and  then  the 
spectrum  as  described  by  Rice  (28). 
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Fig.  ^  ^  Schematic  drawing  of  holocamera:  (A)  dielectric  mirror,  (8)  iris;  (C)  spectrophotometer  cell;  (D) 
ruby  rod-t>4shlamp  assembly;  (E)  iris;  (F)  dielectric  mirror:  (G)  beamsplitter;  (H)  neutral  density  filter;  (I) 
1.6-cm  focal  length  lens;  (J)  25-ji  pinhole;  (K)  43.2-cm  focal  length  collimating  lens;  (L)  front  surface  minor; 

(M)  pin  diode;  (N)  filmhotder  (51) 
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hiqh  intensity  light  such  as  that  available  from  a 
Q-switched.  ruby  laser  is  necessary. 

The  ruby  rod  in  the  laser  is  pumped  by  a  helical  xenon 
flashlamp.  This  assembly  is  located  between  two  flat* 

dielectric  mirrors.  The  optical  duct  formed  has  an  iris  near 
each  mirror;  the  apertures  of  which  can  be  set  to  reduce  the 
number  of  modes  of  the  laser  output.  A  spectroDhotometer 
cell  inside  the  optical  cavity  contains  a  saturable  dye 
called  cr yptocyanine  in  acetone  solvent.  In  a  procedure 
called  O-switchinq.  the  concentration  of  cryptocyanine  is 
increased  allowing  a  reduction  of  the  laser  output  to  a 
single  20  nanosecond  pulse. 

The  laser  output  passes  through  a  beam  expander 

composed  of  a  roughly  1.6  cm  focal  length  lens  which  expands 
the  beam.  a  pinhole  in  its  focal  plane  of  approximately  25 
microns  to  eliminate  i rregul ar i ties  and  a  45  cm  focal  length 
lens  to  collimate  the  expanding  beam.  A  front  surface  mirror 
redirects  the  beam  90  degrees.  The  beam  then  passes  throuqh 
the  water  tunnel  test  section  and  exposes  a  photographic 
plate  held  aqainst  the  test  section  window.  A  neutral 

density  filter  before  the  beam  expander  controls  the  film 

exposure. 

A.2  Reconstruction  System 

For  reconstruction  of  the  test  volume.  the  chemically 
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developed  hologram  is  illuminated  by  another  laser  beam 
through  a  similar  combination  of  expanding  lens*  pinhole* 
and  collimating  lens.  The  holoqram  is  mounted  on  an  XYZ 
vernier  carriage  along  the  axis  of  the  beam.  The 
reconstructed  particles  are  imaged  through  a  microscope  onto 
a  vidicon  and  displayed  on  a  large  close  circuit  video 
screen.  Particles  are  sized  off  the  screen.  Figure  A2  shows 
the  reconstruction  system. 


Figure  A2.  Schematic  Drawing  of  Reconstruction  and 
Viewing  System.  (51) 
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APPENDIX  8 


DATA  ACQUISITION  SYSTEM  SENSITIVITY 
APPLICATION  TO  ACOUSTIC  SPECTRA 

In  order  to  obtain  the  actual  spectrum  of  cavitation 
noi set  a  hydrophone  response  spectrum  had  to  be  considerably 
manipul ated. 

A  receiver  hydrophone  response  spectrum  was  changed 
from  d?V  levels  referenced  to  0*1  volt  to  a  cavitation  noise 
spectrum  by  application  of  a  so  called  sensitivity  spectrum* 
This  spectrum  consisted  of  several  correction  factors  which 
are  included  in  Equation  Bl  to  derive  a  corrective  decibel 
level  for  each  frequency  band* 

dBsens  =  d3v  -  dBcal  -  dBx  ♦  dBd  -  dBm  (Bl) 

where  dBv  is  the  spectral  level  of  the  receiving  hydrophone 
response*  dBcal  is  the  sound  source  strength  per  volt  input* 
d3x  is  the  voltage  input  to  the  sound  source*  dBd  is  a 
correction  for  the  separation  distance  from  the  source  to 
the  receiver*and  dBm  is  a  correction  for  inherent  machine 
attenuati on* 

The  distance  factor  is  qiven  by: 

dBd  =  20  1 oq  (x/  1  meter)  *  dB  re  1  m*  (B2) 
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The  input  voltage  magnitude  correction  is  obtained  from: 


dBx  =  20  log  (V/  1  volt)*  dB  re  1  volt  rms  (B3) 

The  spectral  anal  yzer  had  an  inherent  amount  of  attenuat  i or. 
through  its  electronic  hardware.  This  varied  from  day  to  day 
but  was  on  the  order  of  6  d3V. 

When  all  of  these  factors  are  summed  to  obtain  a 
sensitivity  spectrum*  the  decibel  level  is  in  dB  re  volt  rms 
per  uPa*  In  order  to  obtain  the  correct  final  decibel  levels 
in  dE  re  1  pPa*  some  corrections  had  to  be  made  also  to  the 
spectra  from  the  receivinq  hydrophone  after  spectral 
analysis.  The  corrections  are  given  by  Equation  34. 

dBr  =  dBh  -dBqain  ♦  dBatten  ♦  dBbandwidth  (B4) 

where  dBr  is  the  corrected  received  spectrum  in  dB  re  1  volt 
corrected  by  3  dB  from  peak  to  peak  to  rms*  dBh  is  the 
filtered  receiver  input  to  the  spectrum  analyzer*  dBgain  is 
the  amount  of  gain  used  to  amplify  the  signal*  dBatten  is 
the  amount  of  attenuation  necessary  to  keep  the  spectrum 
analyzer  from  saturating*  dBbandwidth  is  a  correction  of  all 
spectra  to  a  l  Hz  bandwidth.  This  correction  is  calculated 
using  Equation  B5. 

dBbandwidth  =  10  loq  (Af/  1  Hz)  (95) 

where  Af  is  the  frequency  range  for  the  spectrum  divided  by 
512  storage  bins.  While  tnis  is  not  a  totally  correct 
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procedure  because  it  is  oblivious  to  small  frequency  scale 
amplitude  variations*  it  has  been  done  so  that  the  levels 
presented  here  can  be  compared  to  the  levels  obtained  by 
other  investigations  which  also  correct  to  1  Hz  bandwidth 
routinely* 

As  described  in  section  3.6*  a  LC-10  hydrophone  was 
used  as  a  sound  source  next  to  the  model*  The  free  field 
transmitting  response  of  the  LC-10  due  to  a  white  noise 
excitation  was  obtained  in  an  anechoic  calibration  tank  at 
ARL/PSU  usinq  gated  tone  bursts.  This  response  is  shown  in 
Fiqure  51*  Figure  B 2  shows  the  receiving  calibration  for  the 
hydrophone*  .'Jote  that  the  sensitivity  is  fairly  uniform  over 
the  1  -  100  kHz  frequency  range. 

The  response  of  the  receiving  LC-10  hydrophone  to 
acoustic  white  noise  generated  in  the  vicinity  of  the  model 
nose  was  found  by  subtracting  Figure  B1  from  the  received 
spectrum  and  making  the  other  corrections  given  by  Equation 
B 1  to  develop  the  sensitivity  spectrum  of  Fiqure  B3. 

The  receiver  response  to  cavitation  and  tunnel  noise  is 
corrected  with  the  attenuation  and  qain  levels  and  for 
bandwidth*  The  sensitivity  spectrum  was  subtracted  from  this 
receiver  spectrum  to  yield  an  absolute  d9  re  1  pPa  level  for 
the  receiver  response  to  cavitation  noise  spectrum*  Fiqure 
BA  shows  a  typical  final  spectrum. 


115 


Figure  B3. 
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